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Abstract. We highlight recent advances in neutrino astrophysics, the open issues and the 
interplay with neutrino properties. We emphasize the important progress in our understanding 
of neutrino flavor conversion in media. We discuss the case of solar neutrinos, of core-collapse 
supernova neutrinos and of SN1987A, and of the recently discovered ultra-high energy neutrinos 
whose origin is to be determined. 


Neutrinos are weakly interacting particles that travel over long distances and tell us properties 
of the environments that produce them. These elusive particles have kept mysterious for a long 
time. After 1998 many unknown properties have been determined thanks to the discovery of 
neutrino oscillations, first proposed in [I]. This observation was made by the Super-Ivamiokande 
experiment using atmospheric neutrinos [2]. The neutrino oscillation discovery is fundamental 
for particle physics, for astrophysics and for cosmology. 

Neutrino oscillations is an interference phenomenon among the v mass eigenstates, that 
occurs if neutrinos are massive and if the mass (propagation basis) and the flavor (interaction 
basis) do not coincide. The matrix that relates the two basis is called the Maki-Nakagawa- 
Sakata-Pontecorvo matrix [3] which depends on three mixing angles, one Dirac and two 
Majorana CP violating phases (if only three active neutrinos are considered). Solar, reactor 
and accelerator experiments have determined = m 2 — m 2 = 7.6 x 10~ 3 eV 2 , and 

Am 2 2 = m2 — rri\ = 2.4 x 10”' 5 eV 2 , referred to as the atmospheric and the solar mass-squared 
differences respectively. Moreover the sign of Am 2 2 has been measured by the occurrence of the 
Mikheev-Smirnov-Wolfenstein (MSW) effect in the Sun [4j[5j. The sign of Am^ is still unknown 
which makes that there two possible ways of ordering the mass eigenstates. If Am^ > 0 
the lightest mass eigenstate is mi (normal ordering or ’’hierarchy”), whereas if Am 2 ! < 0 
it is m3 (inverted ordering). Most of neutrino oscillation data can be interpreted within the 
framework of three active neutrinos. However a few measurements present anomalies that 
require further clarification. While no single hypothesis of non-standard physics can explain 
all anomalies, the existence of sterile neutrinos, that do not couple to the gauge bosons, can 
partially explain observations (see e.g. [6]). In neutrino oscillation experiments such neutrinos 
manifest themselves through the mixing with the other active species. 

Among the open questions are the mechanism for the neutrino mass, the absolute mass value 
and ordering, the neutrino nature (Dirac versus Majorana), the existence of CP violation in 
the lepton sector and of sterile neutrinos. In the coming decade(s) experiments are planned to 
determine some of these fundamental issues. Cosmological and astrophysical neutrinos span a 
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Figure 1. The left figure shows the proton-proton (pp) nuclear reaction chain with its three 
branches. The pp chain is responsible for energy production in our Sun and low mass stars. The 
theoretical branching percentages define the relative rates of the competing reactions. The right 
figure shows the CN cycle which is thought to play an important role for energy production in 
massive stars. The 15 O and 13 N neutrinos have not been observed yet [13] . 


wide energy range from the meV to the PeV. They are interesting both as probes of the sources 
that produce them and for the information they give us on neutrino themselves (see e.g. HD- 
Historically astrophysical neutrinos have brought milestones in our knowledge of particle physics 
and astrophysics. They can offer avenues for future searches on key open issues. 

1. Solar neutrinos 

Electron neutrinos are constantly produced in our Sun through the proton-proton (pp) nuclear 
reaction chain that produces 99 % of the Sun’s energy by burning hydrogen into helium-4 
j8]. The corresponding solar neutrino flux receives contributions from both fusion reactions 
and beta-decays of 'Be and 8 B. First measured by R. Davis pioneering experiment [9j, solar 
neutrinos were found to be nearly a factor of three below predictions [10] . Over the decades 
solar neutrino experiments have precisely measured electron neutrinos from the different pp 
branches, usually referred to as the pp, pep, 'Be and 8 B and hep neutrinos (Figure 1). The 
measurement of a reduced solar neutrino flux compared to standard solar model predictions - 
the so-called ’’solar neutrino deficit problem” - has been confirmed by these experiments mainly 
sensitive to electron neutrinos, but with some sensitivity to the other flavors. The advocated 
solutions included unknown neutrino properties (e.g. flavor oscillations, a neutrino magnetic 
moment coupling to the solar magnetic fields, neutrino decay, the MSW effect) and questioned 
the standard solar model. 

The MSW effect is a resonant flavour conversion phenomenon due to the neutrino interaction 
with matter. When matter is sufficiently dilute the neutrino interaction with the medium is 
accounted for in the mean-field approximation. The inclusion of the matter potential introduces 
a MSW resonance depending on v properties on one hand (energy, mass-squared difference 
values and signs, mixing angles) and the density profile of the environment on the other. If the 
evolution at resonance is adiabatic, electron neutrinos can efficiently convert into muon and tau 
neutrinos, thus producing a flux deficit. 
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Figure 2. Electron neutrino data on the solar pp, pep, 'Be, 8 B neutrinos from the pp fusion 
chain, as a function of the neutrino energy E v . The results from the Borexino experiment are 
compared to the averaged vacuum oscillation prediction (E u < 2 MeV) and the MSW-LMA 
prediction [E u > 2 MeV), taking into account present uncertainties on mixing angles [18] . 


The solar puzzle was solved by the discovery of the neutrino oscillation phenomenon and the 
results obtained by the SNO and KamLAND experiments (see m for a review). In fact, using 
elastic scattering, charged- and neutral-current neutrino interactions on heavy water, the SNO 
experiment has showed that the measurement of the total 8 B solar neutrino flux is consistent 
with the predictions of the standard solar model : solar electron neutrinos convert into the 
other active flavors. In particular, the muon and tau neutrino components of the solar flux 
have been measured at 5 a mi. Moreover the reactor experiment KamLAND has definitely 
identified the Large Mixing Angle (LMA) solution, by observing reactor electron anti-neutrino 
disappearance at an average distance of 200 km m ■ These observations show that low energy 
solar neutrinos are suppressed by averaged vacuum oscillations while neutrinos having more 
than 2 MeV energy are suppressed because of the MSW effect (Figure 2). Theoretically one 
expects P[y e —> u e ) « 1 — ^ sin 2 26\2 ~ 0.57 (with 0\2 = 34°) for (< 2 MeV) solar neutrinos. 
For the high energy portion of the 8 B spectrum, the matter-dominated survival probability is 
P{v e -> u e ) hi 9 h densit v -> sin 2 0i 2 « 0.31 (see e.g.[33]). 

The Borexino experiment has precisely measured the low energy part of the solar neutrino 
flux, namely the pep pL61, ‘Be m and pp neutrinos JT8]. In fact, by achieving challenging 
reduced backgrounds, the collaboration has just reported the first direct measurement of pp 
neutrinos, the keystone of the fusion process in the Sun. The measured flux is consistent with 
the standard solar model predictions (18). 

The study of solar neutrinos has established that the Sun produces 3.84 x 10 33 ergs/s via the 
pp chain. With the measurement of pp neutrinos, it is clear that energy production in low mass 
stars such as our Sun is well understood. Moreover the occurrence of the MSW effect for high 
energy solar neutrinos shows that different v flavor conversion effects can take place in matter 
compared to vacuum. This phenomenon is also relevant in other contexts, including the early 
universe (at the epoch of the primordial abundances formation), in massive stars, in accretion 
disks around black holes and in the Earth. A future goal for solar experiments is the measurement 









of solar neutrinos from the Carbon-Nitrogen-Oxygen (CNO) cycle which is thought to be the 
main mechanism for energy production in massive stars [S] (Figure 2). Borexino experiment has 
provided the strongest constraint on the CNO cycle that is consistent with standard solar model 
predictions HS3- Moreover precisely determining the transition between the vacuum averaged 
and the MSW-LMA solution is of interest, since deviations from the simplest vacuum-LMA 
transition could point to new physics, such as non-standard interactions [T9], 

2. Supernova neutrinos 

2.1. Core-collapse supernovae and SN1987A 

Supernovae (SNe) are massive stars that undergo gravitational collapse. The implosion of stellar 
cores is at the origin of core-collapse supernovae of type II and Ib/c and of hyper-novae. It was 
realised in [20] that a gravitational binding energy of the order of E « GM^ s /Rns > 10 53 
er£0 associated with the collapse of the star core to a neutron star (NS) would be released as 
neutrino emission and that they eventually play a role in the ejection of the stellar mantle. SNe 
are of type II if they exhibit H lines in their spectra, with type lib having a thin H envelope. 
SNe II-P or II-L present a plateau or a linear decay of the light curves after the peak. SNe are 
of type I when no H lines are seen because the star has lost the H envelope. Type lb shows He 
and Si lines, while Ic shows none of these indicating that the star has lost both the H envelope 
and He shell before collapse. Most of supernova light curves are due to energy deposited by the 
shock or to radioactivity from e.g. ' 56 Ni decay. In some cases the radioactive component can be 
weak or absent. The supernova can still appear as bright if the H envelope is present, otherwise 
it can be invisible (Type Ib/c) (21] . 

The fate of a massive star is mainly determined by the initial mass and composition and the 
history of its mass loss. The explosion produces either neutron stars or black holes, directly or 
by fallback (Figure 3) [21]. SNe initial masses range from 9 to 300 solar masses (M sun ). Stars 
having 9-10 M sun develop an O-Ne-Mg core while if M > 10 M sun an iron core is reached before 
collapse. Hyper-novae are asymmetric stellar explosions with high ejecta velocities, they are very 
bright, producing a large amount of nickel. Hyper-novae are often associated with long-duration 
gamma-ray-bursts. Collapsars are all massive stars whose core collapses to a black hole and that 
have sufficient angular momentum to form a disk (see e.g. mm)- 

On 23 February 1987 Sk -69°202 exploded producing SN1987A, the first naked-eye supernova 
since Kepler’s one in 1604. It was located in the Large Magellanic Cloud, a satellite galaxy of 
the Milky Way. The determined distance is 50 kpc from the Earth based on the expanding 
photosphere method from different groups whose results agree within 10 % (see Table I of [23]). 
This method allows to establish extragalactic distances and cover a wide range, from 50 kpc 
to 200 Mpc. From the observed light-curve and simulations it appears that the core mass of 
SN1987A progenitor was around 6 M sun and total mass ~ 18 M sun and the progenitor radius 
about 10 12 cm [24] . SN1987A is unique because it was observed in all wavelengths from gamma 
rays to radio, and for the first time, neutrinos were observed from the collapse of the stellar core. 
These neutrinos were first discovered by Kamiokande II [25], then by IMB [26] and Baksan m- 
The number of detected electron anti-neutrinos events were 16 in Kamiokande II, 8 in IMB and 
5 in Baksan. Time, energy, SN-angle and background rate for all events are given in the nice 
review [[28]. Several hours before, 5 events were seen in LSD detector that could be due to a 
speculative emission phase preceding the ones seen in the other detectors [29]. Such events are 
often discarded in the analysis of SN1987A data since their are object of debate. The earliest 
observations of optical brightening were recorded 3 hours after neutrinos’ arrival. An enthusiastic 
description of SN1987A discovery is reported in m- 

Three puzzling features concerning SN1987A have set constraints on stellar evolutionary 


1 Mns and Rns are the mass and radius of the newborn neutron star respectively 
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Figure 3. Remnants of massive stars as a function of initial metallicity (y-axis) and initial mass 
(x-axis) [2L. 


models and supernova simulations. The progenitor was a blue supergiant rather than a red 
supergiant, while type II supernovae were thought to be produced by red supergiants. Large- 
mixing processes had transported radioactive nuclei from the deep core far into the H envelope 
of the progenitor and in the pre-supernova ejecta, producing anomalous chemical abundances. 
The presence of three ring-like geometry of the circumstellar nebula around the supernova was 
implying a highly non-spherical structure of the progenitor envelope and its winds (Figure 
[4]) [S3|. Various explanations have been suggested for the presence of these rings, the inner 
one is now dated 20000 years before the explosion. One possibility origin is a binary merger 
event of that epoch M ED- In this case rotation might have played a significant role in 
the dying star. However the prolate deformation of the supernova ejecta at the center of 
the inner ring might have a very different origin (Figure [4]). In fact, the presence of large- 
mixing and of asymmetric ejecta might be suggesting the breaking of spherical symmetry due 
to hydrodynamical instabilities such as the bipolar Standing Accretion Shock Instability (SASI) 
|32j. SN1987A remnant has not been found yet. It is likely not a black hole since the progenitor 
was light enough to be stabilized by nuclear equation-of-states consistent with measured neutron 
star masses [32l 133| . There is currently no sign of a bright pulsar as well, like the one born from 
the supernova explosion in the Crab nebula in 1054. 

Neutrino observations from SN1987A have been used to derive constraints on fundamental 
physics and the properties of neutrinos, axions, majorons, light supersymmetric particles and on 
unparticles. These are derived by the absence of non-standard signatures, by using the intrinsic 
neutrino signal dispersion or the cooling time of the newborn neutron star (see e.g. [Ml). Many 
such limits have been superseded by direct measurements with controlled sources on Earth, 
while others remain. For example, from the three hours delay in the transit time of neutrinos 


























































































Figure 4. Picture of SN1987A 20 years after its explosion with its three rings. The inner 
blowing ring was formed 20000 years before the explosion [36]. 


and photons, a tight limit is deduced on the difference between the neutrino speed c v and light 
speed c, i.e. | (c u — c)/c |< 2 x 10 -9 |35] . 

SN1987A neutrinos have also confirmed the basic features of core-collapse supernova 
predictions concerning the neutrino fluence (time-integrated flux) and spectra. From a 
comparative analysis of the observed neutrino events one gets as a best fit point E = 5 x 1 0 52 
erg and T = 4 MeV for the total gravitational energy radiated in electron anti-neutrinos 
and their temperature respectively (Figure [5J. The emission time is also found to be of 15 
s |28j. According to expectations, 99 % of the supernova gravitational binding energy should be 
converted in v e , v T neutrinos (and anti-neutrinos) in the several tens of MeV energy range. If 
one considers that energy equipartition among the neutrino flavors is rather well satisfied, one 
gets about 3 x 10 53 ergs, in agreement with expectations. Moreover, the neutrino spectra are 
thermal to a fairly good approximation. In this case the average electron anti-neutrino energy 
is related to the temperature by E v = 3 T, giving 12 MeV at the best fit point. For a long 
time this value has appeared much lower than the expected value of about 15 MeV. Currently 
it appears rather compatible with supernova simulations based on realistic neutrino transport. 

Supernova neutrinos are tightly connected with two major questions in astrophysics, namely 
what is the mechanism that makes massive stars explode and what are the sites where the heavy 
elements are formed. Besides the astrophysical conditions and the properties of exotic nuclei, 
neutrinos also determine the heavy element nucleosynthetic abundances. In fact, the interaction 
of electron neutrinos and anti-neutrinos with neutrons and protons in such environments 
determines the neutron-to-proton ratio, a key parameter of the r-process. Among the candidate 
sites for heavy elements nucleosynthesis are core-collapse supernovae, accretion-disks around 
black holes and neutron-star mergers. Numerous studies show that neutrinos impact the neutron 
richness of a given astrophysical environment. In order to finally assess the neutrino impact, 
extensive simulations need still to be performed, which e.g. self-consistently determine neutrino 
and matter evolution, (see the reviews in the Focus Issue [37]). 

Various mechanisms for the SN blast have been investigated in the last decades, including 
thermonuclear, bounce-shock, neutrino-heating, magnetohydrodynamic, acoustic and phase- 
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Figure 5. Energy radiated in electron antineutrinos from SN1987A versus temperature of the 
electron anti-neutrinos. The contours show the allowed regions from a combined analysis of the 
29 events recorded in the three detectors - Kamiokande II, IMB and Baksan. The confidence 
levels are of 1 %, 68.3 %, 90 % and 99 % [281 ). 


transition mechanisms (see [22)). Since the kinetic energy in SN events goes from 10 50 ~ 51 erg for 
SNe up to several 10' 52 erg for hyper-novae, the explosion driving mechanism have, among others, 
to provide such energies. The neutrino-heating mechanism with non-radial hydrodynamical 
instabilities (convective overturn with SASI) appears a good candidate to drive iron-core collapse 
supernova explosions; while the more energetic hypernovae events could be driven by the 
magnetohydrodynamical mechanism |22j. Note that a new neutrino-hydrodynamical instability 
termed LESA (Lepton-number Emission Self-sustained Asymmetry) has been identified [3Sj. 
Successful explosions for two-dimensional simulations with realistic neutrino transport have 
been obtained for several progenitors; while the first three-dimensional explosion has just been 
obtained which shows an enhanced expansion compared to two-dimensional simulations. This 
results from a combination of neutrino heating, SASI and turbulent convection [39] . 

2.2. Neutrino flavor conversion in massive stars 

Important progress has been achieved in understanding how neutrinos change their flavor in 
massive stars. Besides the MSW effect m, calculations have revealed new phenomena due to 
the neutrino-neutrino interaction, the presence of shock-waves and of turbulence (see |4l[fl2] for a 
review). In particular steep changes of the star density profile due to shock-waves induce multiple 
MSW resonances and interference among the matter eigenstates [43* 31]. As a consequence the 
neutrino evolution can become completely non-adiabatic when the shock passes through the 
MSW region. 

The neutrino evolution in flavor space can be understood using the correspondence with 
spin dynamics. In this framework the Schrodinger-like equation that governs the evolution of 
the neutrino flavor amplitudes is replaced by the precession for spins submitted to effective 
magnetic fields (see e.g. [IS]). The components of the magnetic field are the real and imaginary 
part of the off-diagonal term of the Hamiltonian and the difference of the diagonal terms. The 
third component of the spins gives the neutrino flavor content while the other two depend on 
the neutrino mixings. 

Ref. |46 | first pointed out that the inclusion of the neutrino-neutrino interaction introduces 
a non-linear refractive index. Significant effects have been found in gang which has triggered 










intensive investigations since. It has emerged that the neutrino self-interaction produces 
collective stable and unstable modes of the (anti-)neutrino gas. These findings are based on 
the so-called ’’bulb model”, that assumes spherical and azimuthal symmetry of the neutrino 
emission and propagation. Three flavor conversion regimes are found : synchronization, bipolar 
oscillations and spectral split (Figure 6 ). Nearby the neutrino-sphere the neutrino self-interaction 
is so strong that the individual spins are stick together and collectively precess around the 
neutrino self-interaction term. During synchronization no flavor conversion takes place. However 
the corresponding interaction strength decreases with distance r as 1/r 4 [41] . At some point of 
the evolution, the presence of non-zero mixing triggers a flavor instability after which bipolar 
oscillations occur that can be seen either as a pendulum or a gyroscopic pendulum in flavor 
space [49]. In the matter eigenstate basis one can show that the bipolar instability is associated 
with the rapid growth of the matter phase (that depends upon the imaginary components 
of the neutrino self-interaction) [50] . After bipolar oscillations, a spectral split phenomenon 
produces either full or no conversion depending on the neutrino energy. This is a MSW- 
effect in a co-moving frame m, or analogous to a magnetic-resonance phenomenon [52j |. It 
is analogous because the z-axis component of the magnetic field is not constant as is the case 
in the magnetic resonance phenomenon but slowly varies. The collective neutrino modes just 
described produce a swapping of the neutrino fluxes with sharp spectral changes (Figure 7). 
These are flavor conversion mechanisms that emerge in presence of the neutrino self-interaction 
within the ’’bulb model”. The neutrino-neutrino interaction produces interesting conversion 
effects also in accretion disks arising from neutron star-neutron star or black hole-neutron star 
mergers [53] . Note that numerous other oscillation effects have been studied, such as spin-flavour 
oscillations of neutrinos in the generation of pulsar kicks in the gravitation field of magnetized 
neutron stars [53] , 



Figure 6. Electron neutrino survival probability 
in a core-collapse supernova about Is after bounce. 
The neutrino-sphere is taken at a neutron star 
radius of 10 km. The simulation includes 
the neutrino mixing, the neutrino-matter and 
neutrino-neutrino interaction in the single angle 
approximation. The results for a 5 MeV neutrino 
are for inverted hierarchy. Synchronization, 
bipolar oscillations and spectral split can be seen 
in the first 200 km (see text) [ 59] . 



Figure 7. Neutrino fluxes in a core-collapse 
supernova about Is after bounce. The red-dotted 
(black) line corresponds to electron (muon and 
tau) neutrinos at the neutrino-sphere, while the 
yellow (blue dot-dashed) line is the corresponding 
one after the neutrino-neutrino interaction effects. 
The results are at 200 km from the neutrino- 
sphere (at 10 km). The simulation is the same 
as for Figure 6 . Due to the neutrino-neutrino 
interaction, a swapping between the electron and 
the non-electron neutrino types takes place at a 
characteristic energy, here at around 8 MeV [59] . 











Recent investigations have clearly shown that going beyond approximations inherent to the 
’’bulb” model introduces a higher degree of complexity of the flavor evolution. First of all an 
assumption often made within the bulb model is that neutrinos are all emitted at the neutrino- 
sphere with the same angle (’’single-angle” approximation), thus maximizing the coherence of 
the modes. However the interplay between matter and neutrino self-interaction effects need to be 
accurately considered. In fact in the so-called ’’multi-angle simulations” when the matter density 
exceeds or is comparable to the neutrino density, collective effects can be strongly suppressed 
since neutrinos with different emission angles have different flavor histories [55]. In particular it 
appears as for now that collective effects are suppressed during the accretion phase in simulations 
based on realistic density profiles from supernova one-dimensional simulations [56j. If this result 
is maintained in future calculations it implies that during the neutronization burst (50 ms after 
bounce) and the accretion phase (several hundreds ms after bounce) neutrinos change their flavor 
only due to the well understood MSW effect in the outer layers. Note that shock-wave effects 
impact the neutrino fluxes only after about 1 s after bounce in the anti-neutrino (neutrino) sector 
if the neutrino mass ordering is inverted (normal). During the neutronization burst neutrino 
self-interaction effects are absent because, they imply v e v e —> v x v x - In fact, at an early stage of 
the collapse v e only are emitted by electron capture on nuclei. 

Besides the bipolar, other flavor instabilities appear by spontaneous breaking of the 
symmetries of the ” bulb model”, such as the azymuthal symmetry m- This symmetry breaking 
appears even if the initial conditions and the evolution equations have such a symmetry. This has 
been shown using linearised analysis, first proposed in this context in [58] and futher highlighted 
in HP in a broader context. In the small amplitude approximation (also called Random-Phase- 
Approximation in the many-body context) unstable (or runaway) solutions are revealed by the 
presence of complex eigenvalues of the stability matrix. This signals a change in the curvature of 
the energy density of the system around the starting point, given by the initial conditions at the 
neutrino-sphere. Similarly collective stable (or unstable) modes can be identified in many-body 
systems such as atomic nuclei and metallic clusters nn 

The current theoretical description of neutrino propagation in massive stars is based on the 
mean-held approximation. As initial conditions for the neutrino fluxes at the neutrino-sphere, 
results from supernova simulations are taken. The latter is assumed a sharp boundary. For 
example, few collisions outside the neutrino-sphere might modify the picture. Their inclusion in 
a toy model has revealed significant modifications of the flavor patterns m- 

Extended descriptions describing neutrino evolution in dense media have been derived using a 
coherent-state path integral m, the Born-Bogoliubov-Green-Kirkwood-Yvon hierarchy m , or 
the two-particle-irreducible effective action formalism [62] (see also [63 : , B3]). These approaches 
provide a derivation of the evolution equations from first principles and using quantum held 
theorjQ Besides collisions, extended mean-held descriptions present two kinds of corrections : 
spin or helicity coherence [62] and neutrino-antineutrino pairing correlations mi- The former 
are present because of the neutrino mass, while the latter appear in an extended mean-held 
approximation when considering all possible two-point correlators. The most general equations 
for neutrino propagation including both corrections have been derived in m- While the origin 
is very different, both kinds of correlations introduce neutrino-antineutrino mixing, if spatial 
anisotropies of the matter and/or neutrino backgrounds are present. Such corrections are 
expected to be tiny, but the non-linearity of the equations could introduce signihcant changes of 
neutrino evolution in particular in the transition region. This is between the dense region within 
the neutrino-sphere which is Boltzmann treated, to the diluted one outside the neutrino-sphere 
where collective flavor conversion occurs. Numerical calculations are still needed to investigate 
the role of spin coherence or neutrino-antineutrino pairing correlations or of collisions. A first 

2 Note that the quantum field theory treatment of neutrino mixing has been extensively investigated (see e.g. 
|66i [67] I. The possible contribution to the cosmological constant from neutrino mixing is pointed out in [68]. 


calculation in a simplified model shows that helicity coherence might have an impact I55I- 

Another interesting theoretical development is the establishment of connections between 
neutrino flavor conversion in massive stars and the dynamics, or behavior, of many-body 
systems in other domains. Using algebraic methods, Ref. [69] has shown that the neutrino- 
neutrino interaction Hamiltonian can be rewritten as a (reduced) Bardeen-Cooper-Schrieffer 
(BCS) Hamiltonian for superconductivity |70] . As mentioned above, Ref. ni has included 
neutrino-antineutrino correlations of the pairing type which are formally analogous to the BCS 
correlations. The linearisation of the corresponding neutrino evolution equations has highlighted 
the formal link between stable and unstable collective neutrino modes and those in atomic nuclei 
and metallic clusters CD- 

Significant progress has been performed in unravelling how neutrinos change flavor in massive 
stars. However important questions require further investigations, including e.g. the assessment 
of their final impact on a neutrino signal on Earth, of the effect of a wave-packet treatment m 
or of corrections beyond the mean-held approximation EDCDEa, or the implementation of 
realistic matter density profiles in multidimensional simulations and of turbulence. 

2.3. Supernova neutrino observations 

The SuperNova Early Warning System (SNEWS) [73] and numerous other neutrino detectors 
around the world can serve as supernova neutrino observatories if a supernova blows up in the 
Milky Way, or outside our galaxy. Large scale detectors based on different technologies mm 
including liquid argon, water Cherenkov and scintillator are being considered (e.g. JUNO or 
Hyper-K). These have the potential to detect neutrinos from a galactic or an extragalactic 
explosion as well as to discover the diffuse supernova neutrino background produced from 
supernova explosions up to cosmological redshift of 2 (for a review see (76, 77]). 

The observation of the neutrino luminosity curve from a future (extra) galactic supernova 
would closely follow different explosion phases furnishing a crucial test of supernova simulations, 
and information on unknown neutrino properties. In particular, the occurrence of the MSW 
effect in the outer layers of the star and of collective effects depends on the value of the third 
neutrino mixing angle and the neutrino mass ordering. The precise measurement of the last 
mixing angle m EH [80] reduces the number of unknowns. Still, the neutrino signal from a 
future supernova explosion could tell us about the mass ordering, either from the early time 
signal in IceCube ED, or by measuring, in Cherenkov or scintillator detectors, the positron 
time and energy signal associated with the passage of the shock-wave in the MSW region [ 82] . 
Several other properties can impact the neutrino fluxes such as the neutrino magnetic moment 
[83] . non-standard interactions, sterile neutrinos. CP violation effects from the Dirac phase exist 
but appear to be small as established in [82] and further studied in [59] [86] [85]. In spite of the 
range of predictions, the combination of future observations from different detection channels 
can bring key information to this domain (see e.g. [87])■ 

3. Ultra-high energy neutrinos 

The main mission of high-energy neutrino telescopes is to search for galactic and extra-galactic 
sources of high-energy neutrinos to elucidate the source of cosmic rays and the astrophysical 
mechanisms that produce them. These telescopes also investigate neutrino oscillations, dark 
matter and supernova neutrinos (for IceCube). The 37 events collected in IceCube, with 
deposited energies ranging from 30 TeV to 2 PeV, is consistent with the discovery of high 
energy astrophysical neutrinos at 5.7 a (Figure [8]) [SB]. The 2 PeV event is the highest-energy 
neutrino ever observed. 

High-energy neutrino telescopes are currently also providing data on neutrino oscillations 
measuring atmospheric neutrinos, commonly a background for astrophysical neutrino searches. 
Using low energy samples, both ANTARES [89] and IceCube/DeepCore [90] have measured the 
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Figure 8. Arrival directions of the neutrino events seen in the IceCube detectors in Galactic 
coordinates. The grey line shows the equatorial plane. Showerlike events are marked with + 
and muon tracks events with x. A test statistic (TS) for point source clustering at each location 
is shown by the colours. No significant clustering is observed [88] . 


parameters 023 and Ain good agreement with existing data. Neutrino telescope are also 
sensitive to other fundamental properties such as Lorentz and CPT violation m, or sterile 
neutrinos. PINGU, IceCube extension in the 10 GeV energy range, could measure the mass 
hierarchy and be sensitive to the Dirac phase [92] . Such measurements exploit the occurrence of 
the matter effect from neutrinos, both from the MSW and the parametric resonance occurring 
in the Earth [931194] . Feasibility studies are currently ongoing both for PINGU and for ORCA 
[95] to establish if the energy and angular resolution required for the mass hierarchy search can 
be achieved. 
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